Toward Astrophysical Turbulence in the Laboratory 
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Turbulence is a ubiquitous phenomenon in space and astrophysical plasmas, driving a cascade 
of energy from large to small scales and strongly influencing the plasma heating resulting from the 
dissipation of the turbulence. Modern theories of plasma turbulence are based on the fundamen- 
tal concept that the turbulent cascade of energy is caused by the nonlinear interaction between 
counterpropagating Alfven waves, yet this interaction has never been observationally or experimen- 
tally verified. We present here the first experimental measurement in a laboratory plasma of the 
nonlinear interaction between counterpropagating Alfven waves, the fundamental building block of 
astrophysical plasma turbulence. This measurement establishes a firm basis for the application of 
theoretical ideas developed in idealized models to turbulence in realistic space and astrophysical 
plasma systems. 
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Introduction. — Turbulence profoundly affects many 
space and astrophysical plasma environments, playing a 
crucial role in the heating of the solar corona and ac- 
celeration of the solar wind [if, the dynamics of the in- 
terstellar medium 0-0] , the regulation of star formation 
, the transport of heat in galaxy clusters , and the 
transport of mass and energy into the Earth's magneto- 
sphere 0. At the large length scales and low frequen- 
cies characteristic of the turbulence in these systems, the 
turbulent motions are governed by the physics of Alfven 
waves [1] , traveling disturbances of the plasma and mag- 
netic field. Theories of Alfvenic turbulence based on ide- 
alized models, such as incompressible magnetohydrody- 
namics (MHD), suggest that the turbulent cascade of en- 
ergy from large to small scales is driven by the nonlinear 
interaction between counterpropagating Alfven waves 
[l^ . However, the applicability of this key concept in the 
moderately to weakly collisional conditions relevant to 
astrophysical plasmas has not previously been observa- 
tionally or experimentally demonstrated. Verification is 
important because the distinction between the two lead- 
ing theories for strong MHD turbulence ll|, [l^ arises 



from the detailed nature of this nonlinear interaction. 
Furthermore, verification is required to establish the ap- 
plicability of turbulence theories, utilizing simplified fluid 
models such as incompressible MHD, to the weakly col- 
lisional conditions of diffuse astrophysical plasmas. 

Several reasons make it unlikely that the nonlinear in- 
teraction between counterpropagating Alfven waves can 
ever be verified using observations of turbulence in astro- 
physical environments: the spatial resolution achievable 
in astrophysical observations is insufficient, in situ space- 
craft measurements yield information at only a single or 
a few spatial points, and the broad spectrum of turbu- 
lent modes confounds attempts to identify the transfer 
of energy from two nonlinearly interacting Alfven waves 
to a third wave. Only experimental measurements in 



the laboratory can achieve the controlled conditions and 
high spatial resolution necessary. The unique capabilities 
of the Large Plasma Device (LAPD) at UCLA [H], de- 
signed to study fundamental plasma physics processes, 
make possible the first laboratory measurement of the 
nonlinear wave- wave interaction underlying Alfvenic tur- 
bulence. 

This Letter presents the first laboratory measurement 
of the nonlinear interaction between counterpropagating 
Alfven waves, the fundamental building block of astro- 
physical plasma turbulence. The properties of the non- 
linear daugther Alfven wave are predicted from incom- 
pressible MHD theory. The experimental setup and pro- 
cedure are outlined. Analysis of the experimental results 
demonstrate a successful measurement of the nonlinear 
interaction between counterpropagating Alfven waves. 

Theory. — Modern theories of anisotropic Alfvenic 
plasma turbulence are based on several key concepts de- 
rived from the equations of incompressible MHD. These 
equations can be expressed in the symmetric form [3], 



9z±/a<T (vA ■ V)z=' 



-(zT •V)z± - Vp/po, (1) 



where the magnetic field is decomposed into B = Bq + 
(5B, VA ~ 'Bq/ yf^opo is the Alfven velocity due to the 
equilibrium field Bq, p is total pressure (thermal plus 
magnetic), po is mass density, and z^ = v±±S'B±/y/popo 
are the Elsasser fields of the Alfven waves which are in- 
compressible, so that V • z^ = 0. The Elsasser field z+ 
(z~) represents an Alfven wave traveling down (up) the 
mean magnetic field. The second term on the left-hand 
side of ([T]) is the linear term representing the propaga- 
tion of the Elsasser fields along the mean magnetic field at 
the Alfven speed, the first term on the right-hand side is 
the nonlinear term representing the interaction between 
counterpropagating waves, and the second term on the 



right-hand side ensures incompressibility [11 



Consider the nonlinear interaction between two plane 
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Alfven waves with wavevectors ki and k2 , eaeh with non- 
zero components both parallel and perpendicular to the 
equilibrium magnetic field. The mathematical form of 
the nonlinear term in ([1} requires two conditions for non- 
linear interaction to occur: (a) both z+ 7^ and tT ^ 0, 
so the two waves must propagate in opposite directions 
along the magnetic field, implying /cp and k\\2 have op- 
posite signs 0, ; and, (b) the polarizations of the per- 
pendicular wave magnetic fields are not parallel, imply- 
ing k^i X k^2 7^ 0. These properties dictate that the 
fundamental building block of plasma turbulence is the 
nonlinear interaction between perpendicularly polarized, 
counterpropagating Alfven waves. 

For sufficiently small amplitudes, the terms on the 
right-hand side of ([T]) are small compared to the linear 
propa gat ion term, producing a state of weak MHD turbu- 
lence [lO|. Note that the linear term in ([!]) has no coun- 
terpart in incompressible hydrodynamics, eliminating the 
possibility of weak turbulence, a fundamental distinc- 
tion between incompressible hydrodynamic and incom- 
pressible MHD systems. In the weak MHD turbulence 
paradigm, two counterpropagating Alfven waves interact 
nonlinearly to transfer energy to a third wave. This is the 
fundamental interaction underlying the cascade of energy 
to small scales in plasma turbulence. Solving for the non- 
linear evolution using perturbation theory demonstrates 
that the nonlinear three-wave interaction, averaged over 
many wave periods, must satisfy the constraints 
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and oji + u;2 = '^s , 



(2) 



equivalent to the conservation of momentum and energy 



equivale 



17| . Given the dispersion relation for Alfven 
waves, oj = \k^\\vA, and the requirement for counterprop- 
agating waves, the only nontrivial solution to these equa- 
tions has either fcyi = or kp = 0. [l^l 

The frequency constraint in ([2]) ceases to hold when 
the interaction spans only a fraction of a wave period. 
For such brief interactions, energy is transferred nonlin- 
early to a daughter mode at the instantaneous rate given 
by the nonlinear term in ([1]). In the laboratory, this can 
be accomplished by interacting a high-frequency Alfven 
wave ki with a counterpropagating Alfven wave k2 of 
much lower frequency, such that its parallel wavelength 
is much longer than the length over which the waves in- 
teract. The physical effect of the low- frequency wave 
is to generate a shear in the equilibrium magnetic field, 
producing an effective fc||2 = component to the interac- 
tion. The predicted nonlinear product is an Alfven wave 
ka with the properties fcp = fcyi and k^3 = k^i -I- k_L2. 

Experiment. — The first experiment to verify this inter- 
action in the laboratory was performed on the LAPD 
[Tst using a background axial magnetic field of 800 G 
to confine a hydrogen plasma in a cylindrical column of 
16.5 m length and 40 cm diameter. The plasma dis- 
charge exists for approximately 11 ms with a repetition 
rate of 1 Hz. The electron temperature, = 5.0 eV, and 
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FIG. 1. Schematic of the Alfven wave turbulence experiment 
on the LAPD. The Loop antenna generates a large-amplitude 
Alfven wave polarized in the s-direction traveling up the mean 
magnetic field -Bo and the ASW antenna generates a smaller 
amplitude Alfven wave polarized in the y-direction traveling 
down the mean magnetic field. 



density, n = 10^^ cm~'^, were determined using a swept 
Langmuir probe, with a microwave interferometer used 
to calibrate density. The ion temperature, Ti = 1.25 eV, 
was estimated from previous interferometric measure- 
ments. For these parameters, the electron-ion collision 
frequency [l9| is 3 kHz and the characteristic perpen- 
dicular scale is given by the ion sound Larmor radius 
Ps — ^/Te/mi/rii — 0.29 cm, where fli is the ion cy- 
clotron frequency. 

The experimental setup is depicted in Fig. [TJ At one 
end, the Arbitrary Spatial Waveform (ASW) antenna 
"201 21 1 generates an Alfven wave (blue) with frequency 
270 kHz and a wave magnetic field polarized in the 
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y-direction, 6 By, characterized by a precise perpendic- 
ular wavevector ]i±ips = 0.16x. On the other end, a 
Loop antenna [22|, constructed of two crossed current 
loops phased to produce a dominantly horizontal wave 
magnetic field SB^, generates a large-amplitude Alfven 
wave (red) with much lower frequency / = 60 kHz, char- 
acterized by a wavevector dominated by a y-component 
k±2yPs = 0.055. The perpendicular wave magnetic fields, 
SBx and SBy, are measured using an Elsasser probe [23 1 
between the antennas approximately 2 m from the Loop 
antenna and 9 m from the ASW antenna. The measure- 
ments are taken over a perpendicular plane of size 30 cm 
by 30 cm at a spacing of 0.75 cm, yielding a spatial grid 
of 41 X 41 measurement positions. At each position, 10 
shots are averaged to improve the signal-to-noise ratio. 

To generate a nonlinear interaction between the coun- 
terpropagating Loop and ASW Alfven waves, first the 
Loop antenna is turned on and allowed to establish a 
steady wave pattern at / = 60 kHz. Next, the ASW 
antenna launches a perpendicularly polarized, counter- 
propagating wave. The spatial antenna waveform 5B_l, 
measured in single-antenna runs with identical timing, is 
presented in Fig. [5] Colormaps show the magnitude of (a) 
6Bx for the Loop antenna and (b) 5By for the ASW an- 
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FIG. 2. Spatial magnetic-field waveforms from the (a) Loop antenna and (b) ASW antenna. The colormap indicates (a) SB^ 
for the Loop antenna and (b) 5 By for the ASW antenna in mG, and arrows indicate the vector direction of SJi±. 



tenna, and arrows show the vector direction of S'B± . The 
ASW antenna produces a waveform of 30 mG peak-to- 
peak amplitude with little SB^ component, whereas the 
Loop antenna produces a waveform of 3000 mG peak-to- 
peak amplitude that is dominated by the SB^ component. 

The experimental results are presented in Fig. [3] We 
plot (b) the Fourier transform in the perpendicular plane, 
SBx{kxT ky,t), for the Loop antenna by itself and (c) 
5By{kx, ky,t) for the ASW antenna by itself. Note that 
these are the spatial Fourier transforms of the antenna 
patterns in Fig. ^ and both antennas generate a pair of 
waves with ±kj^. 

The analysis of the counterpropagating Alfven wave 
run is designed to identify the nonlinear daughter Alfven 
wave with two distinguishing properties: (a) fc||3 = k\\i, 
which means the daughter wave will have the same fre- 
quency as the ASW wave, and (b) k±^ = k±i + k±2- 
First, to remove the linear contribution to SB^ caused 
by the large-amplitude Loop Alfven wave, we subtract 
the signal of SB^ of the Loop-antcnna-only run from the 
corresponding signal of the counterpropagating run to 
obtain 5Bx(x,y,t). Next, we Fourier transform in time 
the interval with both waves to obtain 5Bx{x, y, /). Since 
the signal of the daughter wave is expected to peak at the 
same frequency / = 270 kHz as the ASW wave, we band- 
pass filter the 5Bx{x, y, f) signal over the frequency range 
170 kHz < / < 370 kHz, setting the Fourier coefficients 
to zero outside this range. The filtered frequency signal 
is then inverse Fourier transformed from frequency back 
to time. Finally, we Fourier transform in the perpendic- 
ular plane to obtain the spatial Fourier transform of the 
daughter wave, SBx{kx,ky,t). 

The result of this procedure, shown in panel (d) of 
Fig. [3l is the key experimental result of this Letter. The 
observational signature of the nonlinear daughter Alfven 
wave is clear, with the wave field SB^ dominated by four 
wavevectors, corresponding to all possible sums of per- 



pendicular antenna wavevectors, k^3 = +k^i ± kj^2 and 
k^3 = — kxi ± kj_2, as shown schematically in panel (a) 
of Fig. 131 The noise level of the subtracted signal in 
Fourier space is 1.2 mG cm^, yielding a signal to noise 
ratio S/N > 10, demonstrating that the measurement 
is physically meaningful. In addition, the amplitude of 
this nonlinear daughter Alfven wave, which agrees to or- 
der of magnitude with theoretical expectations, peaks at 
the ASW wave frequency, /a = 270 kHz, as predicted. 
These results demonstrate that we have successfully mea- 
sured the nonlinear interaction between counterpropa- 
gating Alfven waves, the fundamental building block of 
astrophysical plasma turbulence. 

This experimental finding verifies that the general 
properties of the nonlinear interaction between counter- 
propagating Alfven waves, as derived theoretically in the 
idealized context of incompressible MHD, hold even un- 
der the weakly coUisional conditions relevant to many 
space and astrophysical plasma environments. Although 
these conditions formally require a kinetic description of 
the turbulent dynamics, our results indicate that the key 
concepts derived from reduced fluid models describe the 
essential nature of the turbulent interactions. Future ex- 
periments will probe the turbulent dynamics at smaller 
scales, providing guidance for the extension of existing 
turbulence theories into the uncharted regime of kinetic 
turbulence at scales below the ion sound Larmor radius, 
where the effects of wave dispersion and kinetic wave- 
particle interactions influence the turbulent dynamics. 

Supported by NSF grants ATM 03-17310 and 
PHY-10033446, DOE grant DE-FG02-06ER54890, NSF 
CAREER Award AGS-1054061, and NASA grant 
NNX10AC91G. The experiment presented here was con- 
ducted at the Basic Plasma Science Facility, funded by 
the U.S. Department of Energy and the National Science 
Foundation. 
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FIG. 3. (a) Diagram of kx for the ASW antenna (blue) and Loop antenna (red). For the nonhnear daughter Alfven wave, 
kx3 (black) is the vector sum of the two antenna wave vectors, kxa = +kxi ± kx2 and kx3 = — kxi ± kx2. BuUseyes indicate 
predicted power distribution of the nonlinear product, (b) Colormap of SBx{kx, ky) for the Loop antenna by itself, (c) Colormap 
(mG cm^) of 5By(kx, ky) for the ASW antenna by itself, (d) Colormap of SBx{kx, ky) for the nonlinear daughter Alfven wave. 
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